Realisability of SUSY discrete flavour symmetry in our universe. 
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We study the realisation of supersymmetric discrete flavour symmetry models to the thermal 
history of our universe. We focus on the evolution of the pseudo moduli field among the flavons by 
taking into account finite temperature corrections. We show that the pseudo moduli flavon domi- 
nates the energy density of our universe and this domination makes crucially difficult to realise the 
flavour symmetry models in our universe. We also discuss possible extensions of the supersymmetric 
discrete flavour symmetry models which can ensure the consistency of the models with the thermal 
history of our universe. Finally, we show an extension to realise the thermal inflation by the flavon 
domination. 



1. INTRODUCTION 

Neutrino experimental data provides valuable informa- 
tion on the origin of the observed hierarchies in the lepton 
masses. Recent experiments of the neutrino oscillation 
indicate the almost tri-bimaximal mixing [1] for three 
flavours in the lepton sector, except for large $13 [2]. On 
theoretical aside, many flavour symmetry models have 
been proposed in order to explain the lepton masses and 
mixing angles so far. In particular, non-Abelian discrete 
symmetry, like A4 or S4, have achieved great success in 
deriving the tri-bimaximal mixing [3]. 

In the discrete flavour symmetry models, there are 
a large number of scalar fields, so-called flavon. The 
vacuum expectation values (VEVs) of flavons can pro- 
vide the flavour structures of the standard model (SM) 
fcrmions through the vacuum alignment [4]. It has been 
widely argued that various cosmological issues can be ex- 
plained by the degrees of freedom of flavon fields, for in- 
stance, the leptogenesis scenario with low reheating tem- 
perature [5] or realisation of inflationary universe [6] . 

Especially in the supersymmetric (SUSY) discrete 
flavour symmetry models, there exists a F-flat direction 
among the flavons. It is generally expected that the 
flavour symmetry is spontaneously broken along the flat 
direction among flavons (or equivalently the pseudo mod- 
uli flavon) at sufficiently high temperature in the history 
of our universe. 

In this letter, we study the evolution of the pseudo 
moduli flavon in the SUSY discrete flavour symmetry 
models without relying on specific models. We show that 
the potential energy of the pseudo moduli flavon domi- 
nates the universe at high temperatures, and such domi- 
nation conflicts with successes of the standard Big Bang 
cosmology. Then, we propose a non-trivial extension of 
the models to overcome the problem in the SUSY discrete 
flavour symmetry. With this extension, the domination 
by the flavon potential energy can sufficiently dilute dan- 
gerous relics of the cosmological modulus. Hereafter, we 
denote the pseudo moduli flavon and its VEV as u and 
(u), respectively. 



2. PROBLEM OF PSEUDO MODULI FLAVON 

In a large class of the SUSY discrete flavour symmetry 
models, the pseudo moduli flavon is coupled with right- 
handed neutrinos vr and the VEV of u provides the ma- 
jorana mass term of vrs which is expected to be around 
10 14 GeV to ensure the See-Saw mechanism with order 
one yukawa couplings [7] . We shall consider this class of 
discrete flavour symmetry models below. 

We first give an analytic expression for the VEV of u 
at zero temperature (T = 0) by including SUSY breaking 
effects and 1-loop effective potentials. Next we shall dis- 
cuss the problem of pseudo moduli flavon in the history 
of our universe by taking into account finite temperature 
corrections. 



2.1. Evolution of flavon at T = 

In order to realise the spontaneous symmetry break- 
ing, the flavon potential should have a negative curvature 
around its origin. In the SUSY discrete flavour symme- 
try models at low energies, there is a strong candidate for 
such negative curvature; namely, a tachyonic soft SUSY 
breaking mass of u. 

The 1-loop effective potentials induced by the right- 
handed (s)ncutrinos vr (ur) lift the flavon potential at 
large field value of u. The potential energy of u including 
the soft SUSY breaking mass rh u and 1-loop effective 
potentials at zero temperature can be written by V(u) = 



fh 2 ,\u\ 2 



AV VR + A Vy R , where Vo is the potential 
(u) which is chosen as V((u)) = 0. AV VR 
and AVjjj, are the 1-loop effective potentials from the 
coupling with vr and vr respectively, 
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where M VRl {u) 2 = \y Rl \ 2 \u\ 2 , M 0r1 {u) 2 = \y Rl \ 2 \u\ 2 + 
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rh VR i [8], and the index i represents the generation 
(flavour) of right-handed neutrinos. Here, ym is yukawa 
coupling between u and vri, m 2 VRi is soft SUSY-brcaking 
mass of vru and \i is energy scale where we evaluate 
the potential. 1 We take both driS and m 2 Hi 's to be 
the same order respectively for simplicity. Note that the 
A-terms do not contribute to V(u) because they are for- 
bidden by U(1)r symmetry [4]. We stress that any soft 
SUSY breaking terms do not lift V{u) at large field value 
of u, even if the higher dimensional operators are taken 
into account, in contrast to the analysis in [5]. 2 

By taking the limit of large field value in Eq.(l), i.e., 
\u\ 3> /i, the flavon potential is reduced to 
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Here, we replace the yukawa couplings dri by the av- 
eraged value yn, the validity of which will be discussed 
elsewhere [9]. Consequently, we obtain the analytic ex- 
pression for (u) , 
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: exp 
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We see that the VEV of the pseudo moduli flavon 
strongly depends on the soft SUSY breaking masses and 
the yukawa couplings. We also obtain analytical expres- 
sions for the mass of u around the vacuum and Vq as 
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respectively. 



2.2. Thermal evolution of flavon at T 



Once we take into account the finite temperature ef- 
fect, the quadratic thermal potential T 2 |m| 2 contributes 
to V(u) at high cosmic temperatures. The pseudo mod- 
uli flavon u is trapped at the origin due to this positive 
curvature. Although, one might expect that the discrete 
flavour symmetry is spontaneously broken at sufficiently 
high temperatures, such symmetry breaking is never sim- 
ply realised due to the thermal corrections. 

In the epoch of high cosmic temperatures, the potential 
energy is kept to almost constant until u rolls down to 
its minimum, i.e., 

V(u,T) = V(u,T = 0) + cT 2 \u\ 2 ~V (T > m u ) , (5) 



where the factor c depends on coupling constants in the 
models and is taken to be 0(1) in our following analysis. 

As long as the cosmic temperature satisfies Vq^ 4 > 
T > m u , the flavon potential energy Vq dominates the en- 
ergy density of our universe. We find that the flavon dom- 
ination occurs in many classes of models [10-12], where 
a F-flat direction among flavons exists and the VEV of 
the flat direction is expected to be larger than Aew- We 
will show that achievements of the standard cosmology 
are ruined in consequence of the flavon domination. 

When the cosmic temperature decreases below m u , the 
pseudo moduli flavon u rolls down to the minimum of the 
potential and obtains a non- vanishing VEV, (u) ^ 0. The 
energy of u is simultaneously stored into the SM particles 
through the decays of u. 

In order to recover the standard Big Bang cosmology 
after the flavon domination, the decay of u into the SM 
particles must be completed before the Big Bang nucle- 
osynthesis (BBN). The temperature T R sec ^ of the energy 



transfer is estimated as T R sec ' 

where M p is the reduced Planck mass. We have to build 
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the model to satisfy the condition T R 
to conflict with BBN achievements. 

In order to verify whether y^ sec ) C an be above a few 
MeV, we estimate the decay rate of u. The interactions 
of u in our analysis are given by 



„ ur c r . , Su 
Unt 3 -^Sui/rVr + h.c. + 

V 2 V 2 («) 



Here we decompose u as 
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where a denotes the NG boson associated with the spon- 
taneous breaking of U(l) symmetry in V(u). 3 Using the 
interactions in Eq.(6), the most efficient decay process of 
u into the SM particles is found to be u — > 2v, and the 
decay rate is given by 



T{u^2v) 
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On the other hand, u also decays into the NG boson 
a through the derivative coupling. It is easy to find that 
the decay process u — > 2a dominates the decays of u, and 
the decay rate is given by 



T(u^2a) 



1 
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In our analysis, we will take fi = fh u . 

The author in [">] claimed that the SUSY breaking term, 
A|, w (n 3 +u* 3 )/A, lifts the flavon potential at large field value 
of u. However, such a term is never induced once all of the next 
leading order terms of 1/A is consistently taken into account. 
The relevant terms does not appear due to the definition of F- 
flat direction by including the next leading terms. 



Hence, the most of the energy of u, which dominates the 
energy of our universe, will be stored into the NG boson 



3 In the A4 or S4 symmetry model, this U(l) symmetry can be 
interpreted to the global U(1)b—l symmetry which is implicitly 
included in the supcrpotcntial. 
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a. Once the energy of u is released into a, the energy of 
a is never released into the SM particles before the BBN 
begins, because the interaction rate of a is much smaller 



than one corresponding to T 1 



(sec) 



few MeV. 



According to general discussion given above, we con- 
clude that the SUSY discrete flavour symmetry models, 
for example the models given in [10-12], do not agree 
with the thermal history of our universe, since the en- 
ergy of our universe is overdosed by the NG boson a. 



expect that the model with relatively light right-handed 
neutrinos 4 could be a candidate for the modification of 
the model to overcome the problem of the NG boson a. 
With this extension, the interactions of u are given by 
dnt = ^^x~ ^ uv r v r + h.c. and the branching ratio of 
the process u — > 2a is given by 



Br(u-»2a) 



1 + 24 



1 - 



4m, 
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3. POSSIBLE IMPROVEMENTS 

The essential point of the problem of the overclosure by 
a attributes to the largeness of Br(u — > 2a) and the small- 
ness of the interaction rate of the NG boson a. So long 
as the discrete flavour symmetry is retained rigorously, 
we should extend the models to increase Br(u — > SM), 
since the decay of u into massless NG boson a is always 
opened. 



3.1. Introduction of uQQ coupling 

First we consider the extension by introducing new 
interaction into the superpotential as AW = UquQQ . 
Here, Q and Q arc additional vector-like "heavy quarks" 
which have SU(3)qcd charge. With the aid of this cou- 
pling, the pseudo moduli flavon u can decay into glu- 
ons through the "heavy quark" loops. The decay rate of 
u^rgg is given by [13] 



N 2 a 2 



T (u^-gg) 
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(10) 



where N q denotes the number of flavours of Q, Q. 
The branching ratio of u — > SM is now found to be 



Br(w^SM) 



r(M->gg) 

r(w->2a) 
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To overcome the decay into a, about a hundred "heavy 
quarks" should be introduced. According to the discus- 
sion in [14], however, the running gauge coupling blows 
up due to the corrections from a large number of "heavy 
quarks" below the GUT scale. Thus, this kind of modi- 
fication is not favoured. 



3.2. Light right-handed neutrinos 

Since the process u — > 2v contains the propagators of 
heavy vr, the decay rate of u — > SM is suppressed by 
l/(u). If the pseudo moduli flavon can decay directly 
into the right-handed neutrinos, the branching ratio of 
u — > SM will be enhanced, because the right-handed neu- 
trinos produced by a decay of u subsequently decay into 
active neutrinos and higgses. Thus, one might naively 



According to Eq.(12), Br(tt— > SM) never exceed 50 per- 
cent, even if the right-handed neutrinos are light enough. 
When Br(u — > 2a) is not less than 45 percent, the rela- 
tivistic NG boson a significantly modifies the abundance 
of light elements produced in BBN due to the "speed-up 
effect" [15, 16]. Thus, this kind of modification is not 
efficient for the problem of overclosure by the NG boson. 



3.3. Explicit breaking term of flavour symmetry 

Even if we introduce the new fields which have QCD 
charge, or the sufficiently light right-handed neutrinos, 
the branching ratio of Br(u— >SM) does not mainly con- 
tribute to the total decay width of u. Hence, we consider 
an extension of the model, which kincmatically forbid the 
process u— >2a. In order to give the mass to the NG bo- 
son a, we assume an explicit breaking term for the U(l) 
symmetry in the flavon potential. 5 This also breaks 
discrete flavour symmetry explicitly, as 



AC 



-K (u ■ 



(13) 



where k is a dimensionful parameter. This explicit break- 
ing term gives the NG boson mass m 2 = k/(u). Again, 
we stress that such explicit breaking term cannot be in- 
duced by the SUSY breaking effects or higher dimen- 
sional operators. 

With the help of the explicit breaking term, the de- 
cay process u — > 2a is kinematically forbidden, when k > 
m^(u)/4. To ensure the sizeable T^ ec \ we further in- 
troduce the interaction with vector-like "heavy quarks" , 
AW — VquQQ, with a small number of flavours. 6 The 
energy transfer from the flavon occurs by the process 
u — > gg which can make T R sec ^ around MeV and the width 
of the process is given by Eq.(10). We find that this ex- 
tension is only the possible solution for the problem of 
the overclosure by the NG boson after the flavon domi- 
nation. 



4 The masses of light right-handed neutrinos are induced from the 
higher dimensional operator, i.e., W 9 y ru 2 v c r v c r I 'A [9] 

° Here, we do not specify the origin of the explicit breaking term. 

6 The number of flavour of "quarks" should be less than eight for 
5 representation of SU(5) [14]. 
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4. REALISATION OF THERMAL INFLATION 

As we have already explained, the pseudo moduli 
flavon u is trapped at the origin for T > T cnc i = rh u , and 
the flavon domination occurs when the cosmic tempera- 



ture decreases to T ~ V, 



1/4 



If wc consider the model in- 



cluding the explicit breaking term and the coupling with 
the vector-like "heavy quarks" , the pseudo moduli flavon 
u mainly decays into gluons. It allows us to apply the 
flavon domination to the thermal inflation scenario [17], 
which is the mechanism for diluting the abundance of 
dangerous cosmological modulus X- 

In the presence of thermal inflation, the energy density 
of cosmological modulus is diluted through the entropy 
production after thermal inflation. The modulus abun- 
dance after the thermal inflation is given by [18] 



Px /rpjsec 

s 1 B 



' A 8 R \M P 



(14) 



where A is entropy dilution factor which is defined by 



A 



: J ,(l 



TV(T ond ) T 3 „T<, 



Xo the initial amplitude of cos- 



mological modulus, and Tr the reheating temperature 
after primordial inflation. 

The relevant quantities are estimated as 



-.(sec) 



: exp 
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1 2 Wave 



1 x 10 13 GeV , 



(15) 



90 

7T 2 #» 



1/4 



y/T(u -> gg)M P ~ 2 MeV (16) 

(17) 



where |y avo | and m a 



VR,VQ 



and 



denote the averaged values of 
Here, we set 



m Q ' ™Q j respectively. 
K = m\{u)/A, N q = 8, |j/ ave | = 1, ?7i„ = 500 GeV, and 
Wave = 100 GeV. Eq.(17) indicates that the modulus 
abundance can be sufficiently diluted by thermal infla- 
tion caused by the pseudo moduli flavon u. The details 
of the realisability of thermal inflation by u arc shown 
elsewhere [91. 



5. SUMMARY AND DISCUSSIONS 

We studied the thermal evolution of the pseudo mod- 
uli flavon, which is included in a large number of SUSY 
discrete flavour symmetry models. The pseudo moduli 
flavon has the potential energy due to the finite tem- 
perature corrections, 1-loop effective potentials, and the 
SUSY breaking masses. Especially, we considered the 
model including the coupling between the pseudo moduli 
flavon and the right-handed neutrinos. 

We first gave the analytic expression for the VEV of 
the pseudo moduli flavon at zero temperature. The scale 



of the VEV of the pseudo moduli flavon strongly depends 
on the SUSY breaking masses and yukawa couplings. 

At finite temperature T > rh u , the pseudo moduli 
flavon is trapped at the origin due to the thermal correc- 
tions from the heat bath. During this epoch, the pseudo 
moduli flavon starts to dominate the energy density of 
our universe, when the cosmic temperature decreases be- 
low Vq^ 4 . After the flavon domination ends, the pseudo 
moduli flavon mostly decays into the NG boson asso- 
ciated with the spontaneous discrete flavour symmetry 
breaking. The temperature of energy transfer from a 
into the SM particles is much lower than few MeV, since 
the interaction rate of a is too small. Hence, the energy 
of our universe is overdosed by the NG boson a and the 
realisation of SUSY discrete flavour symmetry models 
becomes difficult significantly. 

To overcome the problem of the overclosure by the NG 
boson, we examined the three possible extensions of the 
SUSY discrete flavour symmetry model : (i) Introducing 
the interaction with vector-like "heavy quarks" into the 
supcrpotential. It opens the decay channel into gluons 
through the "heavy quark" loops. Although this exten- 
sion could throw out the problem with hundred flavours, 
such a large number of "heavy quarks" derive the gauge 
couplings to blow up below GUT scale, (ii) The model 
with light right-handed neutrinos. Although this can 
open the decay channel into the right-handed neutrinos 
of on-shell, the branching fraction of this process cannot 
exceed 50 percent. When the branching ratio of the de- 
cay process, (flavon) — > (SM), is less than 50 percent, the 
relativistic NG boson significantly modify the abundance 
of light nucleus produced in BBN. (iii) Assuming the ex- 
plicit breaking term of the SUSY discrete flavour sym- 
metry in the lagrangian. This breaking term induces the 
NG boson mass and the decay of pseudo moduli flavon 
into the NG boson is forbidden by the sufficiently large 
breaking term. However, the origin of the explicit break- 
ing term is unclear, since such term is never induced by 
any SUSY breaking effects or the higher dimensional op- 
erators, as long as the flavour symmetry is retained ex- 
actly. 

We found that the combination of (i) with the small 
number of flavours of "heavy quarks" and (iii) is the only 
way to avoid the problem of the overclosure by the NG 
boson, and such a combination makes the temperature 
of the energy transfer to be above a few MeV. 

Finally, we showed the applicability of discrete flavour 
model to the thermal inflation scenario. The abundance 
of cosmological modulus is sufficiently diluted by the 
flavon domination. 

We close this letter by pointing out the future tasks on 
our study. In this letter, the values of soft SUSY break- 
ing masses are taken as free parameters. However, if we 
specify the mediation mechanism for the SUSY breaking, 
the size of soft SUSY breaking term would be determined 
by the SUSY breaking scales. The scale of the VEV of 
the pseudo moduli flavon is sensitive to the variance of 
parameters. Therefore the radiative corrections to these 



5 



parameters would significantly change the scale of the 
VEV even if those corrections are small. 

We also comment on another possibility to generate 
the NG boson mass m 2 a in different analysis. Our anal- 
ysis in this letter was based on the F-flat direction and 
unbroken U(1)r symmetry, and that leads to a mass- 
less NG boson a. Once the U(1)r symmetry is bro- 
ken and the nonrenormalisablc Kahlcr potential, such as 
K 3 (u 3 +u^ 3 )/M p for A 4 model and K 3 (u 4 + u^ 4 )/M 2 
for 5 4 model [4, 10, 11], is introduced, the NG boson can 
obtain the mass in analogy with the R-axion [19]. The 
mass of NG boson a in this case is roughly estimated as 
m 2 oc rn 3 / 2 (u) and m 2 oc m 3 / 2 ( u ) 2 /Mp f° r ^4 an d S4 
models respectively. However, within our knowledge, it 
has not been clarified whether the flavour structure of 



the SM fermions is altered by the broken U(l)n symme- 
try and the nonrenormalisablc Kahlcr potential. Detailed 
studies will be given elsewhere [9]. 
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